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Abstract 



Preliminary results from the MILC collaboration for f Ba , fjj t and their ratios 
are presented. We compute in the quenched approximation at /3 = 6.3, 6.0 and 5.7 
with Wilson light quarks and static and Wilson heavy quarks. We attempt to quantify 
all systematic errors other than quenching, and have a first indication of the size of 
quenching errors. 

1 Introduction 

Existing and planned experimental measurements of B-B and B s -B s mixing do not 
constrain the Cabibbo-Kobayashi-Maskawa matrix without knowledge of the heavy- 
light decay constants f B and f Ba and the corresponding 5-parameters. For example, 
the amplitude of B-B mixing is written in the form 

*d ~ \V; d V tb \ 2 f(m t ) ■ % -m 2 B f 2 B B B . (1) 
Here the mixing matrix element is 

\m 2 B f B B B = (5o|(^(l " 7b)6)(^(1 " 7b)*Po) (2) 

and the decay constant f B is defined by the matrix element 

{0\d 7 ^ 5 b\B(p)) =i P J B . (3) 

The only known method of computing f B and the 5-parameter B B from first principles 
is lattice QCD. This fact has led to a major effort in the lattice community to compute 
these quantities [1]. 

The 5-parameter is believed to be close to its vacuum saturation value B B ~ J. and 
the major uncertainty is thus the value of f B . Fortunately this is also the easier matrix 
element to be computed with lattice methods, since it only involves the computation 
of appropriate two-point functions and not the harder three-point functions needed for 
the computation of B B . 

Over the past year and a half, the MILC collaboration has been computing heavy- 
light decay constants in the quenched approximation - the approximation where inter- 
nal fermion loops are ignored - on Intel Paragon computers. Most of the computations 
have been performed on the 512-node Paragon at Oak Ridge National Laboratory, but 
Paragons at Indiana University and at the San Diego Supercomputer Center have also 
been used. More recently we have started computations with the effects of two light 
flavors of dynamical fermions included. This will allow us, eventually, to remove the 
last unknown systematic uncertainty in our computations. 
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2 Features of the lattice computation 



Lattice QCD is a discrete approximation to QCD where the fields are restricted to 
the sites (quark fields) and the links connecting neighboring sites (gauge fields) of a 
space-time lattice with lattice spacing a. Continuum physics is recovered in the limit 
a — > and thus all results have to be carefully extrapolated to this limit. Next, the 
space-time volume is kept finite and the results have thus to be extrapolated to the 
infinite volume limit. Keeping both lattice spacing and volume finite we have a finite 
number of degrees of freedom which allows us to use powerful numerical methods for 
our computations. We use stochastic Monte Carlo methods and therefore our results 
will have statistical errors. 

To calculate we compute, with B an operator that creates a B-meson, e.g., 
B = ^75 6, the two-point function 

_ _ p—E n t 

J2((d 7o75 b)(x,t)B\0)) = £{0|J 7o7B &|n)— -<n|i?t|0) 

2m B V x 1 1 ' y J 

The unknown amplitude (B\B^\0) can be obtained from the two-point function 

(B&) ^ -J—KBlBtlOlV™-'. (5) 

The parameters in eqs. (4,5) are obtained by fits to the (stochastically) computed 
two-point functions. 

Performing such a computation we still do not yet get quite the desired result. 
Though the current we have used for our computation on the lattice looks identical to 
the continuum current, quantum fluctuations induce a finite renormalization 

fF a = Z2 1 fZ, Z A = l + 0(a s ), (6) 

computable in perturbation theory. All our results we will quote include this renor- 
malization. 

For the computation of the two-point functions we need one light and one heavy 
quark propagator. We compute the heavy quark propagator in a "hopping parameter" 
expansion [2] - essentially a 1/mg expansion - keeping 400 terms. This gives very 
good convergence, on our lattices, for quark masses well below the charm. The ability 
to adjust the heavy quark mass arbitrarily is proving very useful in the analysis of 
systematic errors. 

Since we only have results for degenerate light quarks, we determine k s , the strange 



quark hopping parameter, by adjusting the pseudoscalar mass to \J2m 2 K —m 2 , the 
lowest order chiral perturbation theory value. 
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For heavy-light mesons we use the Kronfeld-Mackenzie [3] norm (-\/l — 6k) and adjust 
the measured meson pole mass upward by the difference between the heavy quark 
pole mass ("mi") and the heavy quark dynamical mass ("m 2 ") as calculated in the 
tadpole-improved tree approximation [3]. This procedure eliminates the leading (in 
a s ) discretization errors that would have been of order 0(omg), with amg>l in the 
charm region. 

3 Results 



Table 1: Lattice parameters. 



name 




size 


# configs. 


A 


5.7 


8 3 x 48 


200 


B 


5.7 


16 3 x 48 


100 


E 


5.85 


12 3 x 48 


100 


C 


6.0 


16 3 x 48 


100 


D 


6.3 


24 3 x 80 


100 


F 


5.7; m=.01 
(n F = 2 QCD) 


16 3 x 32 


49 



So far we have results from six different lattices with parameters listed in Table 1. 
j3 = Q/g 2 denotes the bare lattice coupling. 

In the limit of an infinitely heavy quark, heavy quark effective theory tells us that the 
combination (j)p = fp\jMp of pseudoscalar meson decay constant and mass becomes 
independent of the heavy quark mass, up to renormalization group logarithms, with 
corrections that vanish as powers of 1/Mp 

4> P = fp\jM P = (/), 

Therefore, it has become customary to plot fp\jMp vs. 1/Mp. Such a plot for lattice 
D is shown in Fig. 1. The fit is covariant, to the form (7). Although to the eye there 
appears to be reasonable consistency among the heavy-light results and between the 
heavy-light and static-light results, the % 2 /d.o.f for the fit is % 2 (confidence level 
~ 10%), whether or not the static-light point is included. The rather low confidence 
level may perhaps be due to the fact that we have not included additional large-ma 
corrections to the action and operators [4], or simply to the small differences between 
the heavy quark mass and the meson mass Mp. Such effects are under investigation. 
Note that, in an earlier calculation [5], the statistical errors were considerably larger, 



1 + 



M P 



+ 



C 2 

M 2 P 



+ 



(7) 
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Figure 1: fp(Mp)i vs. 1/Mp for lattice D. The light quark is extrapolated to the 
physical mass (m u + m c i)/2. 

and the % 2 /d.o.f for such fits was good. Here the statistical precision has increased to 
a level where small effects are becoming important. 

Table 2 shows results from the six lattices. The lattice-spacing dependence is appar- 

Table 2: Results for decay constants and ratios. /,r=131 MeV scale used throughout. 





A 


B 


E 


C 


D 


F 


Sb 


195(6) 


198(4) 


190(4) 


176(3) 


166(4) 


195(4) 


fBs 


244(5) 


237(3) 


226(3) 


206(3) 


192(3) 


220(4) 


Id 


227(5) 


227(4) 


216(3) 


203(2) 


198(2) 


219(2) 




275(4) 


273(3) 


253(3) 


236(2) 


225(2) 


242(2) 


-f*- 


1.25(2) 


1.20(1) 


1.19(1) 


1.17(1) 


1.16(1) 


1.13(1) 


J p., 

fp 


1.21(1) 


1.20(1) 


1.17(1) 


1.16(1) 


1.13(1) 


1.11(0) 
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ent, but little, if any, finite volume effect is present (compare A and B). This is seen 
more clearly in Fig. 2, which shows /# vs. lattice spacing. It is natural to extrapolate 




0.0 0.5 1.0 

a (Gey) -1 



Figure 2: /# vs. a. Diamonds have scale set by / T ; crosses, by m p . Fits are to the 
diamonds (linear) and to the crosses (linear + quadratic); the a = extrapolations 
are marked by "bursts." The higher cross at a % 0.7 is from lattice B. The "fancy 
diamond" is from configuration F, np = 2 "full QCD," f w scale. 

the quenched /^-scale results linearly to the continuum; we get 148(5) MeV. Note that 
the /^-scale results have much less a dependence than those using m p . This makes 
sense since f w and /# are likely to have rather similar finite a effects. Neither a linear 
nor a linear + quadratic fit to the m p -scale results is particularly good; the linear + 
quadratic fit is somewhat better and gives a result which is consistent, within large 
statistical errors, with the /^-scale results: 158(25) MeV. 

The results on lattice D are consistent with those of [5] which uses the same (3 and 
volume. The major cause of the somewhat smaller central values here is a lower (but 
still consistent) estimate of the scale (1/a ~ 3.0 GeV here vs. % 3.2 GeV in [5]). The 
difference is within the statistical scale error quoted in [5]. 
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We linearly extrapolate to a = all results in Table 2. Systematic errors are then 
estimated — in a preliminary fashion — as follows: 

• Changes of fitting ranges (in t) for the propagators and of types of fits in ljM 
for / pn /Mp give a typical variation of about twice the statistical errors. 

• The dependence on the determination of k s is estimated by finding the change 
in the extrapolated results if k s is fixed using the vector state (j) instead of the 
pseudoscalar state. The difference is especially significant for and is % 0.06 
there. 

• Finite volume effects are estimated by taking the fractional difference between 
results from lattices A and B, using the f w scale. 

• We estimate scale errors by comparing the extrapolated results from the f w and 
m p scales. The difference 10 MeV for the decay constants 0.02 for the 
ratios) is comparable to (but slightly less than) what we would get by comparing 
f w m p -scale values at the smallest available lattice spacing (/3 = 6.3). 

• The effects of using heavy Wilson fermions without the additional corrections to 
the action and operators detailed in [3, 4] are estimated by comparing the original 
fits (see, e.g., Fig. 1) with fits using only the 5 or 6 lightest heavy-light states 
(and, where appropriate, the static-light point). At /3 = 6.3, the maximum value 
of (m 2 — mi)/m 2 is 0.22 with the original fits and 0.04 with the new ones. The 
differences in the results are quite small: % 4 MeV for the decay constants and 

0.01 for the ratios. 

• The effects of quenching are estimated by comparing the results of np = 2 simu- 
lations with the quenched results interpolated to the same lattice spacing. (See 
for instance Fig. 2 for the case of /#.) Note that this seems likely to give an 
overestimate of the quenching effects, since quenched chiral perturbation theory 
arguments [6] indicate that "real" (np = 3) QCD is closer to the quenched the- 
ory than np = 2 QCD. The estimated uncertainty from quenching is of order 
10%, i.e., about 3-4 times the statistical errors, for the decay constants, and of 
order 5% for the ratios of the decay constants. (This is reasonable since some 
cancellations are expected in the ratios.) 

Adding all the above systematic errors except quenching in quadrature, our prelim- 
inary results are 

148(5)(14)(19), f D = 180(4)(18)(16), (8) 
165(4)(20)(17), f Ds = 194(3)(16)(9), (9) 

1.13(2)(9)(4), ^ = 1.09(1)(4)(4), (10) 
Id 



Jb = 

Ib 3 = 

Ib 3 _ 
f B 
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where the decay constants are in MeV. The first error is statistical, the second is due 
to non-quenching systematics, and the third is an estimate of the quenching effect. 
The systematic error estimates on /# and have been recently updated. Further 

study of the systematic errors, and especially the quenching error, is in progress. 

4 Future prospects 

We see that the errors on the heavy decay constants and their ratios are dominated 
by the systematic errors, which are typically a factor of 4-5 larger than the statistical 
errors. Therefore, progress will come from reducing the systematic errors. In particular 
we foresee progress by 

• Removing 0(a) effects by using improved fermionic actions (clover fermions); a 
more ambitious step, correcting for all orders in ma up through order a s , is being 
attempted by the Fermilab group [4] . 

• Diminishing 0(a) effects by using weaker coupling and hence smaller lattice spac- 
ings: work on a simulation with /3 = 6.52 on a 32 3 X 100 lattice has begun. 

• More unquenched results. Since quenching is in many cases the dominant un- 
certainty this appears to be the most fruitful venue. Unfortunately it is also the 
hardest, requiring by far the most computational resources. Existing unquenched 
configurations have relatively large lattice spacings and small volumes, but with 
the advent of more and more powerful computers progress will come. 

We expect that in about two years the systematic errors will be reduced to ±10 MeV 
in fs (not including quenching) with an error of ~ ±12% from quenching. Similarly, 
we expect systematic errors of ±.04 in (not including quenching) plus an error 

of ~ ±.03 from quenching. Since the results seem to depend sensitively on np, the 
prospect of correcting for quenching (as opposed to merely estimating the error) is 
further away: one will need simulations with np = 3 flavors of quarks, two light flavors 
for the u and d quarks, and one somewhat heavier, for the s quark. 

We are grateful to Norman Christ and the Columbia group for providing the np = 2 
QCD configurations, and to Greg Kilcup for help in reading and transferring them. We 
thank A. El-Khadra, A. Kronfeld, and P. Mackenzie for useful conversations. Com- 
puting was done at ORNL Center for Computational Sciences, SDSC, and Indiana 
University. This work was supported in part by the DOE and NSF. 
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